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Method and arrangement for medical X-ray imaging 
BACKGROUND OF THE INVENTION 

5 Three-dimensional X-ray imaging is based on taking several 1-D or 2-D 
. projection images of a 3-D body from different directions. If 1-D projection 
images are available from all around a 2-D slice of the body with dense 
angular sampling, the inner structure of the slice can be determined. This is 
known as Computerized Tomography (CT) imaging technology, which is 

10 widely used in medicine today. A crucial part of CT technology is the 

reconstruction algorithm taking the X-ray images as argument and producing 
a voxel representation of the 3-D body. 

In many practical cases X-ray projection images are available only from a 
15 limited angle of view. A collection of X-ray images of a 3-D body is called 
sparse projection data if (a) the images are taken from a limited angle of 
vieW or (b) there are only a small number of images. Sparse projection data 
does not contain sufficient information to completely describe the 3-D body. 

20 However, some a priori information about the body is typically available 
without X-ray imaging. Combining this information with sparse projection 
data enables more reliable 3-D reconstruction than is possible by using only 
the projection data. 

25 Traditional reconstruction algorithms such as filtered backprojection (FBP), 
Fourier reconstruction (FR) or algebraic reconstruction technique (ART) do 
not give satisfactory reconstructions from sparse projection data. Reasons for 
this include requirement for dense full-angle sampling of data, difficulty to 
use a priori information, for example nonnegativity of the X-ray attenuation 

30 coefficient, and poor robustness against measurement noise. For example 



the FBP method relies on summing up noise elements with fine sampling, 
leading to unnecessarily high radiation dose. 

BRIEF DESCRIPTION OF THE INVENTION 

The aim of the invention is to overcome the problems met in 3-D 
reconstruction of the body that occur when using traditional reconstruction 
algorithms with sparse projection data. This is achieved by a method for 
producing three-dimensional information of an object in medical X-ray 
imaging in which method the object is X-radiated from at least two different 
directions and the said X-radiation is detected to form projection data of the 
object. The object is modelled mathematically independently of X-ray 
imaging. Said projection data and said mathematical modelling of the object 
are utilized in Bayesian inversion based on Bayes' formula 



p(x | m) = 



P pr (x)p(m | x) 



to produce three-dimensional information of the object, the prior distribution 
P pr (x) representing mathematical modelling of the object, x representing the 
object image vector, which comprises values of the X-ray attenuation 
coefficient inside the object, m representing projection data, the likelihood 
distribution p(m|x) representing the X-radiation attenuation model between 
the object image vector x and projection data m, p(m) being a normalization 
constant and the posteriori distribution p(x|m) representing the three- 
dimensional information of the object. 

The invention also relates to a medical X-ray device arrangement for 
producing three-dimensional information of an object in a medical X-ray 
imaging, the medical X-ray device arrangement comprising an X-ray source 
for X-radiating the object from at least two different directions and a 
detector for detecting the X-radiation to form projection data of the object. 
The medical X-ray device arrangement comprises: 



means for modelling the object mathematically independently of X-ray 
imaging 

and means for utilizing said projection data and said mathematical 
modelling of the object in Bayesian inversion based on Bayes' formula 



p(x | m) 



P pr (x)p(m | x) 



p(m) 

to produce three-dimensional information of the object, the prior distribution 
p pr (x) representing mathematical modelling of the object, x representing the 

i 

object image vector, which comprises values of the X-ray attenuation 
coefficient inside the object, m representing projection data, the likelihood 
10 distribution p(m|x) representing the X-radiation attenuation model between 
the object image vector x and projection data m, p(m) being a normalization 
constant and the posteriori distribution p(x|m) representing the three- 
dimensional information of the object. 

15 The invention is based on that with Bayesian inversion it is possible to utilize 
a priori information efficiently in 3-D reconstruction from sparse projection 
data. Any collection of projection data can be used in the reconstruction. 
Application-dependent a priori knowledge can be used to regularize the ill- 
posed reconstruction problem. 

20 

Bayesian inversion improves the quality of 3-D reconstructions over 
traditional methods. In addition, the number of radiographs can be 
minimized without compromising quality of the reconstruction. This is very 
important in medical applications since the X-ray dose of the patient can be 
25 lowered. 

BRIEF DESCRIPTION OF THE FIGURES 

Figurel shows a simple example of X-ray imaging. 

30 
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Figures 2A -2B show a difference between global tomography and local 
tomography. 

Figures 3A - 3B show different types of sparse projection data. Every black 
5 dot represents a location of the X-ray source for taking one 

projection image. 

Figures 4A - 4B show examples of parts of boundaries and cracks of an 
object visible and invisible without a priori information. 

10 

Figure 5 illustrates "pencil beam" X-ray attenuation model. 

Figure 6 is basic flow chart of the method according to the invention. 

15 Figure 7 is an intraoral X-ray device arrangement presenting one 
preferred embodiment of the invention. 

Figure 8 shows a measurement geometry of dental limited-angle 
tomography with a digital intraoral sensor. 



20 



25 



Figure 9 shows results of 3D reconstruction of head phantom in the first 
preferred embodiment. 

DETAILED DESCRIPTION OF THE INVENTION 



In practical imaging situations X-ray images are not always available from all 
around the body. The body might be visible only from certain directions due 
to imaging geometry. For example this is the case in 3-D mammography with 
the breast compressed against the detector, or in intraoral dental imaging 
30 with the detector inside the patient's mouth. This situation is called limited- 
angle tomography. Also, the region of interest might be surrounded by tissue 
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that need not be imaged, like in extraoral dental imaging. This situation is 
called local tomography. In addition the number of radiographs should be 
minimized in medical applications for reducing the X-ray dose of the patient 

5 In the preferred embodiments of the invention Bayesian inversion algorithms 
are used to create a new type of 3-D medical X-ray imaging. It is 
intermediate between a projection radiograph and a full computed 
tomography scan. Two steps are needed to perform successfully such 
imaging: In step one, .the doctor (a) collects all the prior information he or 

10 she has on the tissue of interest, and (b) takes the minimum number of 
radiographs containing the desired 3D information. In step two, a 
tomographic algorithm based on Bayesian inversion is used for 3D 
reconstruction. The algorithm contains two separate mathematical models. 
First, all a priori information (i.e. information that is independent of X-ray 

15 imaging) is used to model the unknown object mathematically. The model is 
put into the form of prior distribution in Bayes' formula. Second, the 
measurement is modelled mathematically. This involves the geometry of the 
imaging device, position of the detector and X-ray source during the 
exposure of each projection image and a model for X-ray attenuation inside 

20 tissue. The mathematical model of the measurement is put into the form of 
likelihood distribution in Bayes' formula. 

In figure 1 is shown a simple example of X-ray imaging, where an X-ray 
source 2 is placed on one side of an object 4 under imaging. Radiation 

25 passes through the object and is detected by a detector 6 on the other side. 
The X-ray source is for example a X-ray source of an intraoral X-ray source 
of a dentist, of a dental panoramic X-ray device, of a surgical C-arm X-ray 
device, of a mammography device or of any other medical X-ray device and 
the detector 6 is a detector of some of those devices. Usually the detector 6 

30 is a digital sensor that can be thought of as a 2-D array almost pointlike 
detectors. 



The 3-D body under imaging is modelled by nonnegative X-ray attenuation 
coefficient. The value gives the relative intensity loss of the X-ray travelling 
within a small distance dr: 



— = —x(r)ar 



5 (1) 

The X-radiation has initial intensity / 0 when entering the object 4 and a 

smaller intensity I x when exiting the object. The following equation shows 
the attenuation law: 

10 f x(r)dr = - [ ^~vdr = log J 0 - log h 

JL Jl I K r ) 

(2) 

where initial intensity /„ is known by calibration, and intensity after object 

I x is known from the corresponding point value in a projection image. Thus 
15 the measured data is the integral of x along the line L. 

In the above model it is not taken into account (a) scattering phenomena 
resulting in X-rays changing their direction, or (b) the dependency of 
attenuation on the X-ray spectrum resulting in low-energy photons being 
20 more easily attenuated than high-energy ones. Effect (b) causes errors in 
measurements and is sometimes referred to as beam hardening. More 
detailed models can be utilized in Bayesian inversion to overcome (a) and 
(b). 



25 



In medical imaging the geometrical arrangements of the X-ray source and 
digital sensor vary according to the diagnostic task and equipment. Figures 
2A-2B and 3 illustrate the types of tomographic data resulting from different 
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imaging situations. For clarity, here are presented two-dimensional 
examples; similar situations can be considered in 3-D. 

In figure 2A-2B are shown two cases according to whether the whole object 
5 4 is fully visible in each projection or not. These cases are called global 

tomography and local tomography, respectively. In figure 2B is radiated only 
ROI (Region Of Interest) 8. 

The types of data described in figures 3A - 3B cover a large range of specific 
10 imaging tasks. The choice of data collection dictates what kind of features 

and details it is easiest to reconstruct reliably from the data. Thus the choice 

of projection views must be made according to the diagnostic task at hand. 

In traditional CT imaging, projections are taken from all around the object 4. 

In the preferred embodiments of the invention radiation dose is lowered by 
15 sampling the angular variable more sparsely. Figure 3A presents more dense 

angular data with full angle and limited angle measurements. Figure 3B 

presents more sparse angular data, also with full and limited angles of 

measurement. 



20 In figures 4A and 4B there are presented examples of parts of boundary and 
cracks that are visible or indetectable in reconstruction without a priori 
information. Figure 4A there is an object 4 under imaging with edge on the 
surface of the object. The edge is detectable on the leftmost position where 
the edge is more parallel to the direction of the X-rays. In the figure 4B there 

25 is a crack inside of the object 4. The crack is detectable in the leftmost 
position parallel to the angle of measurement. 

The main idea in Bayesian inversion method is to consider the inverse 
problem as a problem of statistical inference. All variables are redefined to be 
30 random variables. The randomness reflects uncertainty of their actual values 
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and the degree of uncertainty is coded in the probability distributions of 
these random variables. 

When modelling the measurement mathematically the X-ray attenuation 
5 model and the observation can be assumed linear. The model is presented in 
the following equation. . 

(3) m=Ax+e, 

where the variables m, x and e are random variables. 
10 The matrix A comes from the pencil beam model for the X-ray attenuation. 
This is presented in figure 5. There the unknown 3-D body of the object 4 is 
divided into small voxels 10, and the matrix A contains the lengths of the 
path of the X-ray inside each voxel. In this way the integral in formula (2) is 
approximated with a simple numerical quadrature. 



15 



Assume now that the object image vector x and the noise are jointly 
independent random variables. The conditional probability distribution of x 
given the measurement m is given by Bayes' formula in the form 



p(x m) = — j-^r , 

p(m) 



(4) 

20 where p(m) is a normalization constant. The density p(x|m) is called the 
posteriori distribution of x. p^m-Ax) =p (m|x) is a likelihood distribution 
representing the X-radiation attenuation model between the object image 
vector x and projection data m. 

25 The density P pr , called the prior distribution of x, is designed to contain all 

possible information available of the actual object 4 independently of X-ray 
imaging. It is crucial, in contrast to several classical regularization methods, 
that the choice of the prior distribution should not be based on the projection 
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data. The rule of thumb in the design of prior distributions is that typical 
image vectors (say, of some existing library) should have high prior 
probability while atypical or impossible ones should have low or negligible 
probability. 

5 

In the framework of the Bayesian inversion theory, the posterior distribution 
in formula (4) represents the complete solution of the 3D reconstruction 
problem. 

10 To produce an image of the object 4 based on the posterior distribution, 
several alternatives exist. The most common ones are the maximum a 
posteriori estimator (MAP) and conditional mean estimator (CM). They are 
defined by the formulas 



£>(#map | rn ) = maxjp(a? | ra), 

(5) 

where the maximum on the right hand side is taken over all x, and 



a? CM = J xp{x | m)dx. 



(6) 

Finding the MAP estimator is an optimization problem while finding the CM 
estimator is a problem of integration. 



In figure 6 is presented a basic flow chart of the method according to the 
invention. In method step 600 the object is modelled mathematically 
independently of X-ray imaging. In method step 602 the object is X-radiated 
from at least two different directions. In method step 604 the said X- 
30 radiation is detected to form projection data of the object. In method step 
606 said projection data and said mathematical modelling of the object are 
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utilized in Bayesian inversion based on Bayes formula to produce three- 
dimensional information of the object. The method step 600 is also possible 
to perform after method step 602 or 604. 

In the first preferred embodiment of the invention is presented an application 
to dental radiology. 

X-ray projection images are conventionally used in dental radiology. 
However, certain diagnostic tasks require more precise knowledge of the 3D 
structure of tissue than is available in two-dimensional radiographs. Such 
tasks include implant planning and detection of bone loss between tooth 
roots. 

In figure 7 is presented an intraoral X-ray device 5 arrangement presenting 
the first preferred embodiment of the invention. It is important to notice that 
this is only an example of the medical X-ray device 5 arrangement where the 
invention is possible to be utilized. 

The medical x-ray device 5 in the preferred embodiments of the invention is 
for example a dental panoramic X-ray device, a surgical C-arm X-ray device 
or a mammography device. 

In the first preferred embodiment of the invention the articulated arm 
arrangement 3 moves the X-ray source 2 to the right position. The X- 
radiation begins by pressing the exposure button 12. The X-ray source 2 X- 
radiates the object 4, which is for example teeth of a patient. The detector 6 
detects the X-radiation. The image information which is got by detecting the 
X-radiation is sent by communication link 16 to the computer 14. The 
computer comprises the software means 15 to process the image information 
according to the invention. There can be more than one computer 14 and 
also the software means 15 can situate in more than one computer 14. For 
example the first computer 14 is computer which is used in x-ray imaging. 
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The second computer 14 is computer which is used in processing the image 
information according to the invention. It is possible to have the second 
computer 14 far away from the actual medical x-ray device 5. For simplicity 
in figure 7 is shown only one computer 14. 

In the first preferred embodiment (Fig. 7) of the invention, the dentist's X- 
ray equipment is used for taking a set of 2D projection images that are used 
as input for Bayesian 3D reconstruction algorithm. Such equipment includes 
an intraoral X-ray unit and a digital intraoral sensor. 



Benefits of this approach over conventional CT scan are 

- low cost and convenient usage, 

- high resolution of the projection images, 

- possibility to take as few radiographs as is needed to capture the relevant 
15 3D information, minimizing radiation dose to the patient, 

- the possibility to choose imaging directions so that the X-rays do not pass 
through the whole head but only the interesting tissue, further reducing 
dose. 

20 In the preferred embodiments of the invention, Bayesian inversion is used 
for the 3D reconstruction. Input data for the algorithm is the set of 
projection images and the following a priori knowledge: 
(a) Dental tissue consists of few approximately homogeneous regions with 
well-defined boundaries. 

25 (b) Dental tissue can only attenuate X-radiation, not intensify it. 

The example of the detector 6 used in the first preferred embodiment of the 
invention is based on charge coupled device (CCD) technology and has 
dynamic range of 4096 gray levels. The size of the active imaging area is 
30 34 mm * 26 mm and the resolution is 872 * 664 pixels. After exposure, each 



12 



pixel contains an integer which is proportional to the number of X-ray quanta 
that hit the pixel's area. 

Alternative detectors include any other digital intraoral sensor, digitized X-ray 
film, or any intraoral sensing device converting detected X-ray photons to a 
digital image. 

In the measurement geometry the focal point of the X-ray source 2 moves 
horizontally on a circle in the plane 10, with center at the detector 6, see 
Figure 8. In Figure 8, the detector 6 corresponds to one row of CCD pixels in 
the CCD sensor. 

In figure 8 is presented an example of taking seven radiographs of some 
teeth of a patient. The total opening angle of the projection views is 55 
degrees in this example. Figure 8 represents the xy plane, and the z- 
coordinate axis is thought to be perpendicular to the plane of the paper. 

In the preferred embodiments, mathematical modelling of the object 4 can 
be done with a prior distribution of the form in formula (7): 



where the sum is taken over a collection of 3D neighbourhoods N and the 
value U„(x) depends only on the values of voxels belonging to the 
neighborhood N, and a is a positive regularization parameter used to tune 
the width of the prior distribution. The neighborhood N is typically a set of 
voxels whose center points are closer to each other than some predefined 
maximum distance. 

It is to be noted that formula (7) does not define a probability distribution 
since it does not integrate to 1. Instead, the integral diverges. However, 




(7) 
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when used in the Bayes formula (4) with a suitable likelihood distribution, 
the posterior distribution is integrable. The same remark concerns the priors 
derived from (7), that is> formulae (9), (10) and (11). 

In the preferred embodiments of the invention the three-dimensional 
problem can be reduced to a stack of two-dimensional problems each 
corresponding to a plane determined by a constant value of z. Here figure 8 
represents exactly the situation in the xy plane, i.e. the plane z=0. Each row 
in the detector corresponds to one such 2D problem. This approach leads to 
small approximation error because the X-ray source does not move in the 
correct plane for nonzero z, and this error is neglected. 
Next the modeling of a 2D tomographic problem is explained. Let 



(8) m(j) = A x(j) + e(j) 



denote the jth 2D tomographic problem. Here the vector m(j) contains the 
readings on jth row from each of the seven radiographs. The vector x(j) is 
the jth slice of the 3-D representation x of the object 4 under imaging. The 
matrix A is the pencil beam model matrix for the 2-D tomography problem. 
Now x(j) is a 2D array of pixels. The pixels are denoted by x(j)[k,q], where 
k=l / 2 / ... / K is row index and q=l,2,...,Q is column index. 

The mathematical modelling of the object 4, i.e. incorporation of prior 
information, is next explained for the 2D slice. Define 

K Q 
fc=l q=l 

(9) 

where the value U(N(x(j)[k,q])) depends only on the values of pixels 
belonging to the neighborhood N of the pixel x(j)[k,q] and a is a positive 
regularization parameter used to tune the width of the prior distribution. The 



14 

choice of the neighborhoods N(x(j)[k,q]) is arbitrary, but a typical choice is 
the set of pixels whose centerpoints are closer to the centerpoint of x(j)[k / q] 
than some predefined distance. 

5 In the preferred embodiment, the neighborhoods are chosen to consist of 
two adjacent pixels. Further, the function U is chosen to be a power of 
absolute value of the difference between adjacent pixels. These choices lead 
to formula (10): 
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K-l Q 



A.(» W >- exp(-a(^J2\xW[k,q)-a>W[k + l, q ]\ a + 

fc=l g=l 

+ e e \ xU) m - « w [*. * + 1] n ) 



(10) 

where s is a positive real number and a is a regularization parameter used to 
tune the width of the prior distribution. 

20 In the equation (11) is presented non-normalized total variation (TV) 
distribution, when s=l in equation (10). This is the case according to the 
preferred embodiment of the invention, because when s=l the prior model 
becomes such that it gives high probability for objects consisting of a few 
areas of different attenuation with well-defined boundaries. 



K-X Q 



piTV(*«>) = exp ( - a ( ]T g \x<*>[k, q] - + 1, q] | + 



fc=l .7=1 

30 (11) 



fc=l ?=1 

+ E E h 0) [ fc ' - » w i^9 + 1] I)) 
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X-radiation can only attenuate, not strengthen, inside tissue. This leads to 
positivity prior pos defined by 

(12) pos(x(j)) = 1 if all pixels of x(j) are positive, 0 otherwise. 

5 

The explanation of modeling of one 2D tomographic problem is now 
complete. 

The stack of 2D problems are connected to each other by demanding that 
10 consecutive 2D slices xQ) and x(j-l) should not be very different. 
Mathematically this is expressed as follows: 

(13) pr3D(x(j)) = exp(- yYZ I x(j)[k,q] - x(j-l)[k,q] |), 

15 where the sums are taken over k=l,...,K and q=l,...,Q and y > 0 is another 
regularization parameter. So the prior distribution in formula (4) is the 
product of (11), (12) and (13): 



20 



(14) p pr (x(j)) = prTV(x(j)) pos(x(j)) pr3D(x(j)). 

The measurements are taken into account in the form of the likelihood 
distribution 

p(m& | «W) = cexp{-i(m^ - AjWfir^mW - Ax&)} 
25 (15) 

where £ is the covariance matrix of the Gaussian noise vector e. c ia a 
normalization constant. 



In reality, X-radiation measurement noise is Poisson distributed. Formula 15 
: Y: 30 uses a Gaussian approximation to the Poisson distribution. 
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All the parts are presented of the right hand side of (4). Next it is looked for 
the MAP estimate of x(j), that is, the image x(j) that gives the largest value 
for the posteriori distribution (4). This is equivalent to finding the image x(j) 
that minimizes the expression 

— log(prTV(a?k'>)) - log(pr3D(o;^>)) 

(16) 

with the additional requirement that every pixel of x(j) is positive. 

Minimization of (16) is difficult since F is not a differentiable function due to 
the absolute values present in (11) and (13) and due to the sharp cutting in 
(12). The absolute values in (11) and (13) are replaced by 

15 1*1 ~ M*) = J log(cosh(/?*)), 
(17) 

where 0>O, enabling the use of efficient gradient-based minimization 
routines. The positivity constraint is taken care of by an interior search 
20 minimization method. The results are shown in figure 9 in comparison with 
traditional tomosynthesis technique based on backprojection. 

In this description explained mathematic modelling and other image 
information processing are performed by software means 15 to accomplish 
25 three-dimensional information of the object 4. 

In the second preferred embodiment the invention is utilized in dental 
panoramic x-ray imaging. The object 4 is typically teeth of the patient and 
the medical x-ray device 5 arrangement is dental panoramic x-ray device 
30 arrangement. 
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In the third preferred embodiment the invention is utilized in mammography, 
for example in FFDM (Full Field Digital Mammography). There the object 4 is 
breast of human being and the medical x-ray device 5 arrangement is 
mammography device arrangement. 

5 

In the fourth preferred embodiment the invention is utilized in surgical x-ray 
imaging. There the object 4 is the body part under surgery and the medical 
x-ray device 5 arrangement is surgical x-ray device arrangement. 

i 

10 In the first, second, third and fourth preferred embodiments of the invention 
the basic method steps are same as mentioned with the flow chart in figure 
6. The utilizing of the invention in the second, third and fourth preferred 
embodiment is similar to what is described with the first preferred 
embodiment of the invention and elsewhere in this application except 

15 different medical x-ray imaging applications and their differences because of 
different medical x-ray devices and different objects to be x-ray imaged. 

Although the invention is described above with reference to the examples 
illustrated in the attached figures, it is to be understood that the invention is 
20 not limited thereto but can be varied in many ways within the inventive idea 
disclosed in the attached claims. 
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1. A method for producing three-dimensional information of an object (4) in 
medical X-ray imaging in which method the object is X-radiated from at least 
two different directions and the said X-radiation is detected to form 
projection data of the object (4), characterized in that 

- the object is modelled mathematically independently of X-ray imaging 

- and said projection data and said mathematical modelling of the object are 
utilized in Bayesian inversion based on Bayes' formula 



10 P( x I m ) = 



Ppr (*)Pi™ I x) 



p(m) 

to produce three-dimensional information of the object, the prior distribution 
p P r(x) representing mathematical modelling of the object, x representing the 
object image vector, which comprises values of the X-ray attenuation 
coefficient inside the object, m representing projection data, the likelihood 
15 distribution p(m|x) representing the X-radiation attenuation model between 
the object image vector x and projection data m, p(m) being a normalization 
constant and the posteriori distribution p(x|m) representing the three- 
dimensional information of the object (4). 

20 2. A method according to claim 1, characterized in that the three- 
dimensional information of the object (4) is one or more two-dimensional 
images representing X-ray attenuation coefficient along slices through the 
object. 

25 3. A method according to claim 1, characterized in that the three- 
dimensional information of the object (4) is a three-dimensional voxel 
representation of the X-ray attenuation in the object. 

4. A method according to claim 1, characterized in that the measurement 
30 model is m = Ax+e, where matrix A contains the lengths of the path of the 
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X-ray inside each voxel and the noise e is independent of object image 
vector x leading to the likelihood distribution 

p(m\x) = p noise (m-Ax), 

5. A method according to claim 1, characterized in that the said 
mathematical modelling comprises that X-radiation attenuates when passing 
the object (4), which means that every image voxel is nonnegative. 

6. A method according to claim 1, characterized in that mathematical 
modelling is expressed by the formula: 

P Pr (aO = exp(-a ^ U N (x)) 

N 

where the sum is taken over a collection of 3D neighbourhoods N and the 
value U N _(x) depends only on the values of voxels belonging to the 

neighborhood N, and a is a positive regularization parameter used to tune 
the width of the prior distribution. 

7. A method according to claim 1, characterized in that the 3D 

tomographic problem is divided into a stack of 2D tomographic problems and 
on every such 2D problem, the mathematical modelling is expressed by the 
formula: 

P P r(%) = exp(-a ]$P U N (x)) 

N 

where the sum is taken over a collection of 2D neighbourhoods N and the 
value U N (x) depends only on the values of pixels belonging to the 
neighborhood N, and a is a positive regularization parameter used to tune 
the width of the prior distribution, and the 2D tomographic problems are 
related to each other by the formula 



pr3D(xG)) = exp(- Y ES I x(j)[k /q ] - xG-l)[k /q ] |), 
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where the sums are taken over all pixels (k=l,„.,K, q=l / ... / Q ) and y > 0 is 
another regularization parameter. 

8. A method according to claim 7, characterized in that the neighborhoods 
consist of two adjacent pixels and U calculates a power of the absolute value 
of the difference, leading to the formula 

K-l Q 

P„r(x U) ) = exp ( - o( q]~ + 1, q] \° + 

k=l q=l 

+EEi^ ) [M-* (j ' ) [fc,?+i]n) 

fc=l q=l ' 

where s is a positive real number. 

9. A method according to claim 8, characterized in that s=l corresponding 
to total variation (TV) prior describing objects (4) consisting of different 
regions with well-defined boundaries. 

10. A method according to claim 1, characterized in that mathematical 
modelling is qualitative structural information of the target where the 
structural information is encoded in prior distributions that are concentrated 
around object image vectors x that correspond to the physiological structures 
of the object (4). 

11. A method according to claim 1, characterized in that mathematical 
modelling consists of a list or probability distribution of possible attenuation 
coefficient values in the object (4). 

12. A method according to claim 1, characterized in that the X-ray 
imaging geometry, such as X-ray source position, has unknown error 
modelled in the distribution p(m|x). 
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13. A method according to claim 1, characterized in that the X-radiation 
measurement noise is Poisson distributed. 

14. A method according to claim 1, characterized in that the medical X-ray 
imaging is dental radiography. 

15. A method according to claim 1, characterized in that the medical X-ray 
imaging is surgical C-arm imaging. 

i 

16. A method according to claim 1, characterized in that the medical X-ray 
imaging is mammography. 

17. A method according to claim 1, characterized in that three- 
dimensional information of the object (4) is produced on the basis of the 
maximum a posteriori estimator (MAP) which is calculated by the equation: 

p(#map | rn) = maxp(a? | ra), 

m representing projection data and x representing the object image vector 
and where the maximum on the right hand side of the equation is taken over 



18. A method according to claim 1, characterized in that three- 
dimensional information of the object (4) is produced on the basis of the 
conditional mean estimator (CM), which is calculated by the equation: 



where m represents projection data and x represents the object image 
vector. 



all x. 
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19. A medical X-ray device (5) arrangement for producing three-dimensional 
information of an object (4) in a medical X-ray imaging, the medical X-ray 
device arrangement comprising an X-ray source (2) for X-radiating the object 
from at least two different directions and a detector (6) for detecting the X- 
radiation to form projection data of the object (4), characterized in that 
the medical X-ray device (5) arrangement comprises: 

- means (15) for modelling the object (4) mathematically independently of 
X-ray imaging 

- and means (15) for utilizing said projection data and said mathematical 
modelling of the object in Bayesian inversion based on Bayes' formula 



p(x \ m) 



P pr (x)p(m 1 x) 
p(m) 



to produce three-dimensional information of the object, the prior distribution 
P P r(x) representing mathematical modelling of the object, x representing the 
object image vector, which comprises values of the X-ray attenuation 
coefficient inside the object, m representing projection data, the likelihood 
distribution p(m|x) representing the X-radiation attenuation model between 
the object image vector x and projection data m, p(m) being a normalization 
constant and the posteriori distribution p(x|m) representing the three- 
dimensional information of the object (4). 

20. A medical x-ray device (5) arrangement according to claim 19, 
characterized in that the three-dimensional information of the object (4) is 
one or more two-dimensional images representing X-ray attenuation 
coefficient along slices through the object. 



21. A medical x-ray device (5) arrangement according to claim 19, 
characterized in that the three-dimensional information of the object (4) is 
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a three-dimensional voxel representation of the X-ray attenuation in the 
object. 

22. A medical X-ray device (5) arrangement according to claim 19, 
characterized in that the medical X-ray device arrangement comprises 
means (15) for modelling the measurement as 

m = Ax+e, 

where matrix A contains the lengths of the path of the X-ray inside each 
voxel and the noise e«is independent of object image vector x leading to the 
likelihood distribution 

p(m\x) = p noise (m-Ax). 

23. A medical X-ray device (5) arrangement according to claim 19 
characterized in that the medical X-ray device arrangement comprises 
means (15) for modelling the object (4) mathematically so that X-radiation 
attenuates when passing the object (4), which means that every image voxel 
is nonnegative. 

24. A medical X-ray device (5) arrangement according to claim 19, 
characterized in that the medical X-ray device arrangement comprises 
means (15) for modelling the object (4) mathematically by the formula: 

Ppr(#) = exp(-o; ^2U N (x)) 

N 

where the sum is taken over a collection of 3D neighbourhoods N and the 
value U N (x) depends only on the values of voxels belonging to the 
neighborhood N, and a is a positive regularization parameter used to tune 
the width of the prior distribution. 
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25. A medical x-ray device (5) arrangement according to claim 19, 
characterized in that the 3D tomographic problem is divided into a stack 
of 2D tomographic problems and on every such 2D problem, and the medical 
X-ray device arrangement comprises means (15) for modelling the object (4) 
mathematically by the formula: 

P»r(x) = exp(-at^TU N (x)) 

N 

where the sum is taken over a collection of 2D neighbourhoods N and the 
value U N (x) depends only on the values of pixels belonging to the 

neighborhood N, and a is a positive regularization parameter used to tune 
the width of the prior distribution, and the 2D tomographic problems are 
related to each other by the formula 

pr3D(x(j)) = exp(- y££ | x(j)[k,q] - x(j-l)[k,q] |), 

where the sums are taken over all pixels (k=l,...,K, q=l,...,Q ) and y > 0 is 
another regularization parameter. 

26. A medical X-ray device (5) arrangement according to claim 25, 
characterized in that the neighborhood systems consist of two 
neighboring pixels xj, xk or voxels xj, xk and U N (x) calculates a power of the 

K-l Q 

p pr (* 0) ) - exp ( - a( \*<*>[k t q] - + l,q)\ 3 + 

fc=l q=l 
fc=l 9=1 

absolute value of the difference, leading to the formula where s is a positive 
real number and a is a regularization parameter used to tune the width of 
the prior distribution. 
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27. A medical X-ray device (5) arrangement according to claim 26, 
characterized in that the medical X-ray device arrangement comprises 
means (15) for modelling the object (4) mathematically by setting s=l 
corresponding to total variation (TV) prior describing objects consisting of 
different regions with well-defined boundaries. 

28. A medical X-ray device (5) arrangement according to claim 19, 
characterized in that the medical X-ray device arrangement comprises 
means (15) for modelling the object (4) mathematically by assuming that 
mathematical modelling is qualitative structural information of the target 
where the structural information is encoded in prior distributions that are 
concentrated around image vectors x that correspond to the physiological 
structures of the target. 

29. A medical X-ray device (5) arrangement according to claim 19, 
characterized in that the medical X-ray device arrangement comprises 
means (15) for modelling the object (4) mathematically by assuming that 
mathematical modelling consists of a list of possible attenuation coefficient 
values in the object. 

30. A medical X-ray device (5) arrangement according to claim 19, 
characterized in that the medical X-ray device arrangement comprises 
means (15) for modelling the object (4) mathematically by assuming that the 
X-ray imaging geometry, such as X-ray source position, has unknown error 
modelled in the distribution p(m|x). 

31. A medical X-ray device (5) arrangement according to claim 19, 
characterized in that the medical X-ray device arrangement comprises 
means (15) for modelling the object (4) mathematically by assuming that X- 
radiation measurement noise is Poisson distributed. 
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32. A medical X-ray device (5) arrangement according to claim 19, 
characterized in that the medical X-ray imaging is dental radiography. 

33. A medical X-ray device (5) arrangement according to claim 19, 
characterized in that the medical X-ray imaging is surgical C-arm imaging. 

34. A medical X-ray device (5) arrangement according to claim 19, 
characterized in that the medical X-ray imaging is mammography. 

35. A medical X-ray device (5) arrangement according to claim 19, 
characterized in tiiat the medical X-ray device arrangement comprises 
means (15) for producing three-dimensional information of the object (4) on 
the basis of the maximum a posteriori estimator (MAP), which is calculated 
by the equation: 

p(^MAP | m) = maxp(# | ra), 

m representing projection data and x representing the object image vector 
and where the maximum on the right hand side of the equation is taken over 
all x. 

36. A medical X-ray device (5) arrangement according to claim 19, 
characterized in that the medical X-ray device arrangement comprises 
means (15) for producing three-dimensional information of the object (4) on 
the basis of the conditional mean estimator (CM), which is calculated by the 
equation 



where m represents projection data and x represents the object image 
vector. 




ABSTRACT 

The invention relates to a medical X-ray device 5 
arrangement for producing three-dimensional 
information of an object 4 in a medical X-ray imaging 
medical X-ray device arrangement comprising an X-ray 
source 2 for X-radiating the object from different 
directions and a detector 6 for detecting the X-radiation 
to form projection data of the object 4. The medical X- 
ray device 5 arrangement comprises: 

- means 15 for modelling the object 4 mathematically 
independently of X-ray imaging 

- and means 15 for utilizing said projection data and 
said mathematical modelling of the object in 
Bayesian inversion based on Bayes' formula 



p(x\m) = 



P pr ^)p{m\x) 



p{m) 

to produce three-dimensional Information of the object, 
the prior distribution p pr (x) representing mathematical 
modelling of the object, the object image vector x, 
which comprise values of the X-ray attenuation 
coefficient inside the object, m representing projection 
data, the likelihood distribution p(m|x) representing the 
X-radiation attenuation model between the object 
image vector x and projection data m, p(m) being a 
normalization constant and the posteriori distribution 
p(x|m) representing the three-dimensional information 
of the object 4. 
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the object is modelled mathematically 
independently of X-ray imaging 



\/ 



the object is X-radiated 

from at least two different directions 
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^ the said X-ray radiation is detected 
604~\ to form projection data of the object 



\/ 



said projection data and said mathematical modelling 
of the object are utilized in Bayesian inversion based on 
Bayes formula to produce three-dimensional information 
of the object 
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